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1. Introduction 



Gauge fields can be used to model the long distance behaviour of several condensed 
matter systems 0, a connection which has been particularly exploited for planar systems 
. In a nutshell, the idea is that charge fluctuations around a given ground state are 
decribed by a conserved current j^, which in (2+1) dimensions can be represented in 
terms of a gauge field according to oc e^°"^dotB^. For a wide class of systems 
the effective action governing the dynamics of the charge fluctuations is quadratic in the 
gauge fields at long distances [|l|. Clearly this effective action is also gauge invariant, 
reflecting the original gauge invariance of the definition of the current j^: one obtains 
thus an effective gauge theory at long distances (which is not necessarily relativistic). 
The ground states of a wide class of planar condensed matter systems @ can thus be 
classified according to the lowest derivative term appearing in their effective gauge theory 
at long distances. This way Chern-Simons terms describe incompressible quantum fluids 
(quantum Hall states) and chiral spin liquids while the Maxwell term describes a (2- 
dim.) superfluid (superconductor) 0. 

In this paper we shall investigate a further connection between Abelian gauge theories 
and certain condensed matter systems, namely Josephson junction arrays . In a recent 
publication [H we studied non-perturbative features of the (2+l)-dimensional gauge theory 
with mixed Chern-Simons term 

C = -^F.^F'^' + ^A,e^-^d^B^ - ^U.F^ (1-1) 

when the gauge symmetries associated with the two Abelian gauge fields and B^ are 
compact. We also pointed out the relevance of (|1.1| ) to the zero-temperature physics of 



planar Josephson junction arrays. Here we will derive and study this connection in detail. 

After reviewing in section 2 the basic physics of ( p. . 1| ) and our lattice notation, we 
shall show in section 3 that the zero-temperature partition function of Josephson junction 
arrays in the self-dual approximation coincides with the Euclidean partition function of 



the lattice version of ( p. . 1|) with periodic mixed Chern-Simons coupling. This means that 
the two gauge fields are compact variables only as far as their coupling is concerned. The 
periodicity is implemented by two types of topological excitations which constitute 
electric and magnetic closed loops with short-range interactions. The two energy scales 
of Josephson junction arrays, the charging energy Ec and the Josephson coupling Ej are 
directly related to the two massive parameters and of ( |1 . 1| ) . 
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In section 4 we investigate the non-perturbative structure of this Chern- Simons lattice 
gauge model. The phase structure is determined by the 3-dimensional statistical mechan- 
ics of the topological excitations and reflects the self-duality of the model. We find three 
possible phases at zero temperature. For small e/ g there is a superconducting phase with 
logarithmic confinement of magnetic fluxes; in this phase the original Ra gauge symmetry 
of (|1 . 1|) is broken down to Za so that the full symmetry is given by Za x Rb- Corre- 
spondingly, one of the two massive excitations of ( |1 . Ij ) becomes massless. The dual phase 
is realized for large e/g. In this phase we have logarithmic confinement of electric charges 
and symmetry Ra x Zb, with a corresponding massless excitation. An infinite energy (volt- 
age) is required to separate the charge dipoles and produce a current through the sample: 
we call this phase with infinite resistance a "superinsulator". Depending on the details 
of the lattice, a third phase can open up between the superconductor and the superinsu- 
lator. In this phase the topological excitations are irrelevant, the symmetry is Ra x Rb 
and both excitations are massive. The amount of energy required to produce a current 
through the sample is exponentially small. In |^ we called this phase the Chern-Simons 
phase: in presence of dissipation it would actually correspond to a "metallic" phase of the 
model . The superconductor-insulator quantum phase transition is actually observed 
experimentally in planar Josephson junction arrays at very low temperatures . 

Recently it has been suggested that Josephson junction arrays in presence of Uq offset 
charges and external magnetic fluxes per plaquette might have quantum Hall phases 
for either charges fl^ or fluxes [|12| , depending on the ratios rig/n^ and Ec/Ej. In 
section 5 we discuss these purely two-dimensional quantum Hall states in the framework 
of the gauge theory representation. Specifically, we show that they can be described by 
additional pure Chern-Simons terms for either one of the two gauge fields or S^. In 
this phases the charges and vortices combine to form an incompressible quantum fluid 
of charge-fiux composites with short-range interactions. Localized excitations are charge 
and fiux carrying anyons [^| . 

We then investigate how one of the distinctive feature of Josephson junction arrays, 
namely the periodicity of charge-vortex couplings affects these quantum Hall states. We 
find that this periodicty can induce two types of phase transitions. The charge-fiux fiuid 
corresponding to the charge quantum Hall phase can either expel the fiux and form a charge 
superfiuid corresponding to a conventional superconductor or condense into a charge-fiux 
superfiuid. Correspondingly, the fiux-charge fiuid corresponding to the vortex quantum 
Hall phase can either expel the charge and form a fiux superfiuid corresponding to a 
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superinsulator or condense into a flux-charge superfluid. These superfluids of charge-flux 
composites are (logarithmic ) oblique confinement phases [|T3[ []T^ corresponding to anyon 
superconductors fl^ . We thus conclude that Josephson junction arrays might provide the 



first explicit realization of the anyon superconductivity mechanism. 

In section 6 we generalize our results to three dimensions (even if three-dimensional 
Josephson junction arrays have not yet been fabricated). In this case, one of the two 



gauge fields becomes an antisymmetric Kalb-Ramond tensor gauge field [|T8[ and the (3-M) 



dimensional gauge theory we obtain is the so-called BF-model . This is an Abelian 
gauge model with a conventional Maxwell gauge field and a Kalb-Ramond gauge field cou- 
pled by a topological mass term. In three dimensions the magnetic topological excitations 
become compact surfaces on the lattice and self-duality is lost. The zero-temperature phase 
structure is determined by the statistical mechanics of a model of coupled random loops 
and random surfaces in four Euclidean dimensions: this can also be viewed as the Euclidean 
partition function for a lattice model of particles interacting with closed Nielsen-Olesen 
type strings. While the statistical mechanics of random loops is by now well developed 
1^ there is no corresponding amount of analytical results for random surfaces 

. Assuming three distinct phases as in (2+1) dimensions, with condensation of electric 
loops, no condensation of topological excitations and condensation of magnetic surfaces 
we can identify the first two again with superconducting and metallic phases, respectively. 
In the phase with condensation of magnetic surfaces the charge dipoles are bound by 1/r 
potentials, which are long-range but not confining. Therefore only a finite amount of en- 
ergy is required to separate them and the system behaves as an insulator (as opposed to 
a superinsulator in two dimensions). 

2. The lattice Chern-Simons model 

Our model ( p. . 1|) can be rewritten in terms of the dual field strengths 



1 

r = 

as follows 



[2.1) 



Throughout this paper we use units such that c = 1 and h = 1. 
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For later convenience we have introduced a magnetic permeability ry, equal for the two 
gauge fields. The coupling constants and g'^ have dimension mass, whereas the coefficient 
K of the mixed Chern-Simons term is dimensionless. Note that we take to represent a 
pseudovector gauge field, so that the mixed Chern-Simons term does not break the discrete 
symmetries of parity and time reversal. 

The action corresponding to (p.2|) is separately invariant under the two Abelian gauge 
transformations 



(2.3) 



with gauge groups Ra and Rb, respectively. Moreover, the action is also invariant under 
the duality transformation 

A^^B^ , 

(2.4) 

so that the model is self- dual. 

The Lagrangian ( |2.2|) can be easily diagonalized by the linear transformation 

(2.5) 

B^ = J- {a^ - bfj,) . 

In terms of these new variables the model (|2.2D describes a free theory, 

Ccs = -- (-GoG^ + G,G') + ^a^e^«^a«a, - - (-gog^ + g^g') - ^h^e^^^d^h, , 
eg \r] ) 27r eg \ri J 2tx 

(2.6) 

where G^ and g^ are the dual field strengths for the new gauge fields and h^j,, respectively. 



This Lagrangian describes a doublet of excitations with topological mass [52 



m = ^ , (2.7) 



and spectrum 



£;(q) = . /to2 + -|q|2 . (2.^ 



In the following we shall formulate a Euclidean lattice version of the above Chern- 
Simons model. To this end we introduce a three-dimensional rectangular lattice with 
lattice spacings l^ in the three directions. In particular we shall take the lattice spacings 



li = I2 = I and identify /q witli tlie spacing in tlie Euclidean time direction. Lattice 
sites are denoted by the three-dimensional vector x; the gauge fields A^{x) and B^{x) are 
associated with the links (x, fx) between the sites x and x + ft, where ft, denotes a unit 
vector in direction /i on the lattice. 

On the lattice we introduce the following forward and backward derivatives and shift 
operators: 

d,f{x) ^ f(^ + ^>^f')-f(^) , s,f{x) ^ fix + l,f) , 

(2.9) 



f{x) - f{x - l^f) 

I 



d^f{x) = S^f{x) = f{x-l^fi). 



Summation by parts on the lattice interchanges both the two derivatives (with a minus 
sign) and the two shift operators; gauge transformations are defined using the forward 
lattice derivative. Corresponding to the two derivatives in ( |2.9| ), we can define also two 
lattice analogues of the Chern-Simons operators e^auda 0: 

kfj,u = S nd^Qii/doi , k^i/ = dfictudctSi/ , (2.10) 

where no summation is implied over equal indices /i and u. Summation by parts on the 
lattice interchanges also these two operators (without an extra minus sign) . The operators 
( pj.lO| ) are both local and gauge invariant, in the sense that 

kfiudi/ = d^k^i/ = , k^jydi/ = d^k^i/ = , (2.11) 

and their product reproduces the relativistic, Euclidean lattice Maxwell operator: 

k^akav ~ knakctiy — ^/xj^V + d^d^ , (2.12) 

where = d^d^ is the three-dimensional Laplace operator. Using k^^, we can also define 
the lattice dual field strengths as 

F — h A 

(2.13) 

f = k Ft 

The identity ( |2.12| ) then tells us that we can simply write the relativistic, Euclidean lattice 
Maxwell terms as Yjx -^a^-^m ^'^d Yjx ft^ffj-- 
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Using all these definitions we can now write the Euclidean lattice partition function 
of our model (p.2|) as follows: 



^2.14) 



Zcs = j'^^f^J ^M-Scs) , 

where we have introduced the notation VAf^ = Y[{x fi) dAfj,{x) and gauge fixing is under- 
stood. 

For later convenience we introduce also the finite difference operators 

= lf,d^ , Af, = l^d^ , (2.15) 

where no summation over equal indices is implied. Correspondingly, we introduce also the 
finite difference analogue of the operators k^^ and k^^j\ 

Kfj^i, = S'^c^q-i/Aq, , K = e^Q,;y Ao.iS'jy . (2.16) 

These satisfy equations analogous to ( p.ll|) and ( |2.12|) with all derivatives substituted by 
finite differences. 



3. Josephson junction arrays 

Josephson junction arrays H] are quadratic, planar arrays of spacing / of superconduct- 
ing islands with nearest neighbours Josephson couplings of strength Ej. Each island has 
a capacitance Cq to the ground; moreover there are also nearest neighbours capacitances 
C. The Hamiltonian characterizing such systems is thus given by 

X <xy> ^ ^ 

where boldface characters denote the sites of the two-dimensional array, < xy > indicates 
nearest neighbours, is the electric potential of the island at x and $x the phase of its 
order parameter. For generality we allow for any integer in the Josephson coupling, so 
that the phase has periodicity 2n/N: obviously N = 2 for the real systems. 
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With the notation introduced in the previous section the Hamiltonian (|3.1j ) can be 
rewritten as 



H^J2 lv{Co-CA)V + J2 Ej (1 - cos N (A,$)) , (3.2) 

where A = A^A^ is the two-dimensional finite difference Laplacian and we have omitted 
the exphcit location indeces on the variables V and $. 

The phases $x are quantum-mechanically conjugated to the charges Qx on the islands: 
these are quantized in integer multiples of (Cooper pairs for N = 2): 

Q = QeNpo , po e Z , (3.3) 



where is the electron charge. The Hamiltonian ( |3.2| ) can be expressed in terms of charges 
and phases by noting that the electric potentials are determined by the charges Qx via 
a discrete version of Poisson's equation: 

(Co - CA) Fx = Qx . (3.4) 

Using this in ( p^) we get 

H^J2 ^^^c Pog^^Po + 5^ (1 - COS iV(A,$)) , (3.5) 

X C x,i 

where Ec = q1/2C. The integer charges po interact via a two-dimensional Yukawa po- 
tential of mass ^JC^jC jl. In the nearest-neighbours capacitance limit C ^ Co, which is 
accessible experimentally, this becomes essentially a two-dimensional Coulomb law. From 
now on we shall consider the limiting case Cq = 0. In this case the charging energy Ec 
and the Josephson coupling Ej are the two relevant energy scales in the problem. These 



two massive parameters can also be traded for one massive parameter y/2N^EcEj, which 
represents the Josephson plasma frequency and one massless parameter Ej /Ec- 

The zero-temperature partition function of the Josephson junction array admits a 
(phase-space) path-integral representation . Since the variables po are integers, the 



imaginary-time integration has to be performed stepwise; we introduce therefore a lattice 
spacing Iq also in the imaginary-time direction. This has to be just smaller of the typical 



time scale on which the integers po vary, in the present case the inverse of the Josephson 
plasma frequency: lo <0 (l / \/2N'^EcEj^ . We thus get the following partition function: 

J—k/N 

(3.6) 



{Po} 



S = J2-^N PoAo$ + N^Eclo Po^Po + Yl ^oEj (1 - cos N (A,$)) , 

X x,i 

where now the sum in the action 5" extends over the three-dimensional lattice with spacing 
Iq in the imaginary time direction and / in the spatial directions. 

In the next step we introduce vortex degrees of freedom by replacing the Josephson 
term by its Villain form : 

/+Tr/N 
exp(-5) , 

{-a (3.7) 
S = J2-iN poAo$ + N^Eclo Po^Po + ( + ^^^] ■ 

Strictly speaking, this substitution is valid only for IqEj ^ 1; however the Villain approx- 
imation retains all most relevant features of the Josephson coupling for the whole range of 
values of the coupling Ej and therefore we shall henceforth adopt it. 

We now represent the Villain term as a Gaussian integral over real variables pi and 
we transform also pq to a real variable by introducing new integers vq via the Poisson 
summation formula 

k=-\-oo n=+oo 

^ eicp{i2Tikz) = ^ S{z - n) . (3.8) 

fc= — oo n= — oo 

By grouping together the real and integer p and v variables into three- vectors p^ and Vfj,, 
// = 0, 1, 2 we can write the partition function as 

/■ p+tt/N 

Z=y2 VpJ eM-S) , 

r J-Tt/N 

S = Y, -^Np, (a,$ + ^z;^) + N'Eclo Po^Po + ^ • 

Following we use the longitudinal part of the integer vector field to shift the 
integration domain of $. To this end we decompose as follows: 

= A^m -f A^a + K^^ip^, , (3.10) 
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where m E Z,\a\ <1 and K^^ defined in ( p.l6|) . Here the vectors ^ip^ are not integer, but 



they are nonetheless restricted by the fact that the combinations = K^,^Vi, = K^aKavi^u 
must be integers. The original sum over the three independent integers {f^} can thus be 
traded for a sum over the four integers {m, g^} subject to the constraint A^g^ = 0. 
The sum over the integers {m} can then be used to shift the $ integration domain from 
[— tt/A^, +7r/A^) to (—00, +00). The integration over $ is now trivial and enforces the 
constraint A^p^ = 0: 



(3.11) 



{<?m} 

S = y^ -«27r p^K^.i;, + NHqEc Po —rPo + 



We now solve the two constraints by introducing a real gauge field 6^ and an integer 
gauge field a^: 

p^ = K^i^hj, , e R , 

(3.12) 

Qfi = Kf,^a„ , a^e Z . 

By inserting the first of these two equations and by summing by parts, the first term in the 
action ( |3.11| ) reduces to Ylx ~^27r6^q'^. By inserting the second of the above equations and 
by summing by parts again, this term of the action finally reduces to the mixed Chern- 
Simons coupling — z27r a^K^^h^. Using the Poisson formula ( P7^ ) we can finally make 
also real at the expense of introducing a set of integer link variables {Q^} satisfying 
the constraint A^Q^, which guarantees gauge invariance: 



{Qm} 



^ 1 p^ ^'-''^ 



In this representation K^^b^, represents the conserved three-current of charges, while K^ud^ 
represents the conserved three-current of vortices. Note that, actually, both these con- 
served currents are integers (the factors of are explicit): indeed, the summation over 
{Q^} makes (and therefore also K^^^a^) an integer, and then the summation over {a^} 
makes K^jjb^, an integer. The third term in the action (|3.13|) contains two parts: the 



2 

longitudinal part (pf ) describes the Josephson currents and represents a kinetic term for 



the charges; the transverse part [pf) can be rewritten as a Coulomb interaction term for 
the vortex density qo by solving the Gauss law enforced by the Lagrange multiplier 60 • 

The partition function ( p.l3| ) displays a high degree of symmetry between the charge 
and the vortex degrees of freedom. The only term which breaks this symmetry (apart from 
the integers Q^) is encoded in the kinetic term for the charges (Josephson currents). This 
near-duality between charges and vortices has already been often invoked in the literature 
1^ to explain the experimental quantum phase diagram at very low temperatures. Here 
we introduce what we call the self-dual approximation of Josephson junction arrays. This 
consists in adding to the action in ( |3.13| ) a bare kinetic term for the vortices * and combining 
this with the Coulomb term for the charges into J2x nHqEc '^^^ coefficient is chosen so 
that the transverse part of this term reproduces exactly the Coulomb term for the charges 
upon solving the Gauss law enforced by the Lagrange multiplier oq. The longitudinal part, 
instead, represents the additional bare kinetic term for the vortices. Given that now the 
gauge field has acquired a kinetic term, we are also forced to introduce new integers 
via the Poisson formula to guarantee that the charge current K^j^h^, remains an integer: 



ZsD = ^ "^^M / exp(-S's£,) 



{Q^} 

2 22 



z z z 

SsD = V -^27r a^K^^K + 777%^ + i^rrh- + ^STra^Q^ + i2T:h^M^ 

^ ZIqJ^j iV^io-C/C 



where the new integers satisfy the constraint A^M^ = to guarantee gauge invariance. 
After a rescaling 

^0 = —?=—ao , = ^ , 

VNlo VNl 
o - 27r , _ 27r , ^^.15) 

we obtain finally 

ZsD= Yl J ^^A* / ^^/^ ^M-Ssd) , 



lol^ ^ ^ Jo/'« . , „ y . . (3-16) 



+ (loQoAo + IQ^Ai) + i^ (loMoBo + IM,B, 



* Note that such a kinetic term is anyhow induced by integrating out the charge degrees of 
freedom. 
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where and fi are defined in ( |2.13| ) and 

e''=2NEc, K = N, g^ = ^Ej. (3.17) 

This is exactly the partition function of our lattice Chern-Simons model ( p.l4|) in the limit 
of infinite magnetic permeability 77 = oo and with additional, integer-valued link variables 
and coupled to the two gauge fields. Note that, with the above identifications, the 
topological Chern-Simons mass (|2.7|) coincides with the Josephson plasma frequency: 



m 



^2N^EcEj . (3.18) 



In the physical case N = 2 this reduces to m = ^/8EcEJ. From the kinetic terms in ( |3.16| ) 
we can also read off the charge and vortex masses: 

_ 1 _ N 

"^*=Pg2 = 2NPEc ' 

In the regime ml < 0(1), which is typically experimentally relevant, we can choose 
Iq = I: in this case the infinite magnetic permeability constitutes the only non-relativistic 
effect in the physics of Josephson junction arrays in the self-dual approximation. However, 
we expect this non-relativistic effect to be irrelevant as far as the phase structure and the 
charge-vorticity assignements are concerned. Therefore, for simplicity, we shall henceforth 
consider the relativistic model, by setting Iq = I and i] = 1, although it is not hard to 
incorporate a generic value of 77 into our subsequent formalism: 



ZsD= Yl j ^^m/ eM~SsD) 



{Qm} 

{M^} (320) 



SsD = ^ -^F^F^ - i'^-^A^k^^B^ -f- -^f^f^ + U^/kA^Q^ + ily/KB^^M^ . 

Josephson junction arrays in the self-dual approximation constitute thus a further, exper- 
imentally accessible example of the ideas presented in and 0. The action in ( |3.2q ) 
provides in fact a pure gauge theory representation of a model of interacting charges and 
vortices, represented by the conserved currents 

3 

charge _ ^ r. rj 

27r (321) 

/vortex _ ^ t A 
11 



^ ^ J-'k , , „ l' 



where the prefactors are chosen so that the quantum of charge is given by k, while the 
quantum of vorticity is given by 1/k (factors of and 27r are absorbed in the definitions 
of the gauge fields and the coupling constants). 

In this framework, the mixed Chern-Simons term represents both the Lorentz force 
caused by vortices on charges (coupling of q'^^^^s'^ to the "electric" gauge field A^) and, 
by a summation by parts, the Magnus force ||2^ caused by charges on vortices (coupling 



of (pjl^^^^^ to the "magnetic" gauge field -B^). The integer- valued link variables and 
represent the (Euclidean) topological excitations in the model. They satisfy the 
constraints 

d^^Q^^ = , 



d^M^ = . 



(3.22) 



In a dilute phase they constitute closed electric (Q^) and magnetic (M^) loops on the 
lattice; in a dense phase there is the additional possibility of infinitely long strings. Due 
to the constraints ( |3.22| ) we can choose to represent these topological excitations as 

Q jj. = Ikfii/Yiy , Yj^ G Z , 



(3.23) 



and reabsorb them in the mixed Chern-Simons term as follows: 

^ E ■ ■ ■ - i^^..) -^y.)^-- (3.24) 

In this representation it is clear that the topological excitations render the charge-vortex 
coupling periodic under the shifts 

2% 



2. (3.26) 

In physical terms, the topological excitations implement the well-known periodicity of 
the charge dynamics under the addition of an integer multiple of the fiux quantum 1/k 
per plaquette and the (less-known) periodicity of the vortex dynamics under the addition 
of an integer multiple of the charge quantum k per site. 

If we would require that the full action (including charge-charge and vortex-vortex 
interactions) ( |3.2C1| ) be periodic under the shifts ( |3.25| ), then we would obtain the compact 
Chern-Simons model studied in 0. In this case the relevant topological excitations would 
be essentially iX^ and lY^: since these can also describe finite open strings, there is the 
additional possibility of electric and magnetic monopoles [|[. As we showed in [^, these 
monopoles play a crucial role in the regime ml ^ 1. 
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4. Phase structure analysis 

In this section we investigate symmetry aspects and non-perturbative features of the 
model ( |3.2C1| ) due to the periodicity of the charge-vortex interactions encoded in the mixed 
Chern-Simons term. As expected, these depend entirely on the topological excitations 
which enforce the periodicity. 

Upon a Gaussian integration the partition function (|3.20| ) factorizes readily as 

ZsD = Zcs ■ Ztop , (4.1) 
where Zcs is the pure gauge part defined in ( p.l4|) and 



Ztop = ^ exp (-5'top) , 

{Qm} 

Stop ^ E ^ Q^J^^-^ + t ''^J^''^ 

X 

-I- ? O ^ M 

^ I '^'^V2(m2-V2) 

with m defined in ( |2.7| ), describes the contribution due to the topological excitations. The 
phase structure of our model is thus determined by the statistical mechanics of a coupled 
gas of closed or infinitely long electric and magnetic strings with short-range Yukawa 
interactions. The scale (1/m) represents the width of these strings. In our case it is of the 
order of the lattice spacing /. The third term in the action ( |4.2|) , describing the topological 
Aharonov-Bohm interaction of electric and magnetic strings, vanishes for strings separated 
by distances much bigger than (1/m): in this case the denominator reduces to m^V^ and, 
by using either one of the two equations in ( p.23| ) and the constraints ( |3.22| ) one recognizes 



immediately that the whole term in the action reduces to (i27rinteger) , which is equivalent 
to *. 

4.1. Free energy arguments 

In order to establish the phase diagram of our model we use the free energy arguments 



for strings introduced in [^^l ^^^d extensively used for the analysis of four-dimensional self- 



dual models 16 



* This reflects the fact that the original charges and vortices satisfy the Dirac quantization 
condition. 
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The usual argument for strings with Coulomb interactions |27] is that interactions 
between strings are unimportant for the phase structure because small strings interact via 
short-range dipole interactions, while large strings have most of their multipole moments 
canceled by fluctuations. This argument is even stronger in our case, where the interaction 
is anyway short-range. Therefore one retains only the self-energy of strings, which is 
proportional to their length, and phase transitions from dilute to dense phases appear 
when the entropy of large strings, also proportional to their length, overwhelms the self- 
energy. We shall also neglect the interaction term between electric and magnetic strings 
(imaginary term in the action ( [4.2|) ). This is clearly a good approximation if both types 
of topological excitations are dilute. 

Thus, one assigns a free energy 

^^(m/) Q2 + ^^Giml)M^ -^N (4.3) 



to a string of length L = IN carrying electric and magnetic quantum numbers Q and M, 
respectively. Here G{ml) is the diagonal element of the lattice kernel G{x — y) representing 
the inverse of the operator P (m^ — V^). Clearly G(ml) is a function of the dimensionless 
parameter ml. The last term in ( [4.3|) represents the entropy of the string: the parameter 
11 is given roughly by /i = ln5, since at each step the string can choose between 5 different 
directions. In ( [4.3|) we have neglected all subdominant functions of A^, like a InA^ correction 
to the entropy. 

The condition for condensation of topological excitations is obtained by minimizing 
the free energy ( [4.3| ) as a function of A^. If the coefficient of in ([4.3D is positive, the 
minimum is obtained for A^ = and topological excitations are suppressed. If, instead, the 
same coefficient is negative, the minimum is obtained for N = oo and the system will favour 
the formation of large closed loops and infinitely long strings. Topological excitations with 
quantum numbers Q and M condense therefore if 

le^KGiml) ^9 Iq'^KGiml) , ,9 , , . 
1 — L q2 ^ Ji V — ; < 1 . (4.4) 

If two or more condensations are allowed by this condition one has to choose the one with 
the lowest free energy. 

The condition ( |4.4| ) describes the interior of an ellipse with semi-axes 2ij,/{le'^i^G{ml)) 
and 2n/{lg'^KG{ml)) on a square lattice of integer electric and magnetic charges. The 
phase diagram is obtained by investigating which points of the integer lattice lie inside the 
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ellipse as its semi-axes are varied. We find it convenient to present the results in terms of 
the dimensionless parameters Im and e/g: 

mlG{ml)TX J I < 1 , electric condensation , 

H ^ ^ 1 I > 1 , magnetic condensation , 

- < — ,^ , electric condensation , /a rc\ 

mlG{ml)n ^ I u ^ e ^ mlG{ml)-K 1 ,. 

^ ' > 1 — > < — ,N < - < — — , no condensation , 



e ^ mlG(ml)TT , • i , • 

- > — — , magnetic condensation 



As expected, these condensation patterns are symmetric around the the point e/g = 1, 
refiecting the self-duality of the model. In first approximation the electric (magnetic) 
condensation phase is characterized by the fact that {Q^} ({M^}) fiuctuate freely, while all 
= {Q^ = 0). Within this approximation it is clearly consistent to neglect altogether 
the interaction term between electric and magnetic strings in ( |4.3|) . Taking into account 
small loop corrections [|2^ in the various phases can lead to a renormalization of coupling 
constants and masses and, correspondingly, to a shift of the critical couplings (mZ)crit and 
{e/g)crit for the phase transitions. A notable exception is the case in which there is only 
one phase transition: in this case the critical coupling is {e/g)cvit = 1 due to self-duality. 

4-2. Wilson and 't Hoop, loops 

In order to distinguish the various phases we introduce the typical order parameters 
of lattice gauge theories ||2^ , namely the Wilson loop for an electric charge q and the 
Hooft loop for a vortex cj): 



Lvi/ = exp z-^ Iq^A^ ) 



(4.6) 

Lh = exp I {(pK^ ^ /(/>^S^ 

where and vanish everywhere but on the links of the closed loops, where they take 
the value 1. Since the loops are closed they satisfy 

df^q^^ = di^(l)f, = . (4.7) 

The expectation values {Lw) and {Lh) can be used to characterize the various phases. 
First of all they measure the interaction potential between static, external test charges q 
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and —q and vortices (j) and — ^, respectively [|[. Secondly, by representing the closed loops 
and (p^ as 

q^ = Ik^jyAI , 



(4. 

Ik 



we can rewrite the Wilson and 't Hooft loops as 



(4.9) 



which is a lattice version of Stoke's theorem, the integers A'^ and A'j^ (= 0, ±1) representing 
the area elements of the surfaces spanned by the closed loops. The second terms of the 
expansions of (Lw) and (Lh) in powers of q and (j) measure therefore the gauge invariant 
correlation functions {Ffj_{x)F^{y)) and {f^{x)fi,{y)). Third, if we represent (j)^ as 

<f><f>^ ^ ^ k^^Al-^- , (4.10) 
we can also rewrite the 't Hooft loop as 

Lh = exp ^^ . (4.11) 

With the interpretation of A^™' as an external electromagnetic gauge potential the expec- 
tation value of the 't Hooft loop measures the electromagnetic response of the system in 
the various phases. An analogous relation clearly holds for the Wilson loop. 

The expectation values of the Wilson and 't Hooft loops are easily obtained by com- 
bining the definitions ( [4.6|) with (|3.20|) : 



^ Ztop(Q^,M^) 



(4.12) 

(Lh) 



Ztop(Q^,M^) 

where the notation is self-explanatory. In the following we shall analyze these expressions 
in the various phases obtained in ( [4.5|) . We shall mostly only indicate the form of small 
loop corrections: a full renormalization group analysis is beyond the scope of the present 
paper and we won't be able to predict the orders of the phase transitions. 
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Let us begin with the electric condensation phase. In this phase the ground state 
contains many infinitely long electric strings Q^. These have a crucial effect on the gauge 
symmetry associated with the gauge field A^. To see this let us consider a gauge transfor- 
mation Afj^ —* + d^A, where, for simplicity, we take A as a function of the component 
only. If we choose the usual boundary conditions = = at infinity, the change 
of the action ( p.20|) under the above gauge transformation is given by 

ASsD = ^iV^ (A(x^ = +oo)Qi(x^ = +00) - A(x^ = -oo)Qi(x^ = -00)) . (4.13) 

In a dilute phase, with only small closed loops, Qi(x^ = +00) = Qi{x^ = —00) = 
and the action is automatically gauge invariant. In a dense phase, with many infinitely 
long strings, Qi{x^ = +00) and Qi{x^ = —00) are generically different from zero. Gauge 
invariance requires that ASsd vanishes modulo i2n. In the dense phase this is realized 
only if A takes the values 

271- 

A = —n , n e Z , (4.14) 

^/^^ 

at infinity. This means that, in the electric condensation phase, the global gauge symmetry 
is spontaneously broken down to the discrete gauge group Z, so that the total (global) 
symmetry of this phase is Za x Rb- 

The Wilson loop expectation value takes a particularly simple form if the external test 
charges are multiples of the charge quantum: q = tik, n E Z. In fact, since we sum over 
{Qm}' integer loop variables ng^ can be absorbed by a redefinition of the appropriate 
Q^'s, with the result 

{Lwiq = nK)) = 1 . (4.15) 

This indicates that, in this phase, external test charges q = nn are perfectly screened 
by the topological excitations and behave thus freely. In order to compute the Wilson 
loop expectation value for generic q we have to perform explicitly the sum over {Q^}- 
To this end we have to remember the constraint d^Q^ = 0. We solve this constraint by 
representing = Ik^^nv and summing over {n^}, with the appropriate gauge fixing. 
We then use Poisson's formula ( |3.8| ) to turn this sum into an integral, by introducing a 
new set of integer link variables {k^} satisfying d^k^ = in order to guarantee the gauge 
invariance under + Id^i. At this point we can perform explicitly the Gaussian 

integration over {n^j^}. In the approximation of neglecting terms proportional to V^/m^ 
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(keeping such terms would not alter substantially the result) the new integers {kfj_} can be 
absorbed by a redefinition of the magnetic topological excitations {M^}, giving the result: 

/ r \ Zcorr. (?/i) 

[Lw) ~ 



■^corr. (Qfi — 0) 

^corr.(g,)= e^PE(-§^/^Z^^-+^27r^A^M^V 

{Mp} loops X ^ ^ 

Since the magnetic topological excitations are in a dilute phase we have to sum only over 
small closed loops: in this phase the dominant part of ln{Lw) vanishes for generic q and 
the whole result is given by small loop corrections. These are identical in form to the 
small loop corrections for the correlation functions in the low-temperature phase of the 



three-dimensional XY model [25|; correspondingly the Wilson loop expectation value can 



be computed by exactly the same low-temperature expansion used for the XY model [2^ . 



The first-order term in this expansion is obtained by considering only the smallest possible 
lattice loops and gives the result 

{Lw) = exp I 2e-^ [cos (s^r^c/^) " l] ) • (4-17) 

The periodicity of this result is a direct consequence of the spontaneous symmetry breaking 
Ra Za- This implies also that the gauge invariant correlation function reduces to 

{F^{x)FM) « ('^A^.V - rf^rf.) -jf , (4.18) 

which is essentially a contact term on the scale of the lattice spacing. 

The computation of the 't Hooft loop expectation value follows exactly the same lines 
as the above computation of the Wilson loop. The results is 



{Lh) = exp -— 2^ 




(4.19) 



The first few terms in the expansion of the small loop corrections can again be computed 
with the same techniques as in the low-temperature phase of the XY model ||25|| . One 
finds that their contribution amounts to perturbative corrections of the Coulomb coupling 
constant g^K^cf)^ /2l of the dominant term in ( |4.19| ). 
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From ( [4.19| ) we can extract the nature of the electric condensation phase. First of aU, 
by considering, as usual, a rectangular loop of length T in the imaginary time direction and 
of length R in one of the spatial directions and computing the dominant large-T behaviour 
of \n{LH) we find that the interaction potential between external test vortices of strength cj) 
and —(f) is proportional to livR. Vortices are thus logarithmically confined, which amounts 
to the Meissner effect. Secondly, by using the representations ( [4 .81) and ( ^^91) , we find the 
correlation function 

{U^)My)) oc ^ , (4.20) 

which is long-range, indicating that the "S^-photon" is niassless. This is the massless ex- 
citation associated with the spontaneous symmetry breaking of the global gauge symmetry 
Ra Za- Third, by using the representations (|4.10|) and (|4.11|) , we find that the induced 
electromagnetic current is given by 

Jr ^ - ^] , (4.21) 

which is the standard London form. We thus conclude that the electric condensation phase 
is actually a superconducting phase. 

No further computation is needed to extract the nature of the magnetic condensation 
phase: this is the exact dual of the electric condensation phase just described. Specifically, 
the global gauge symmetry associated with is spontaneoulsy broken down to Zb, so 
that the total symmetry of this phase is Ra x Zb. Correspondingly, the "A^-photon" is 
massless and the {F^{x)F,y{y)) correlation function is long-range. Electric charges are log- 
arithmically confined, which means that an infinite energy (voltage) is required to separate 
a neutral pair of charges. We call this phase with infinite resistance a superinsulator. In 
real Josephson junction arrays we expect however the conduction gap to be large but finite 
due to the small ground capacity Co, resulting in a normal insulator. 

If mlG{ml)Tx / li > 1 a third phase can open up between the superconducting and 
superinsulating phases. In this third phase both the electric and the magnetic topological 
excitations are dilute. Far away from the phase transitions and to first approximation we 
can neglect them altogether. This gives the result 



IT \ ( ^^^^ ^i^^ \ 

{Lh) = exp 



21 ^""m^ 
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(4.22) 



Small loop corrections to these results can be obtained by restricting the {Qfj,} and {M^} 
sums in ( |4.12| ) to small closed loops and using again the same techniques as in the low- 
temperature expansion of the XY model These will lead to perturbative corrections of 
the coupling constants and masses in (|4.22|) ; however the first-order result ([4.22| ) is enough 
to establish the nature of this phase. The global symmetry characterizing this phase is 
Ra X Rb and, corrrespondingly, both "photons" are massive, resulting in short-range 
correlation functions {F^{x)F^j(ii)) and {f^{x)fv{y)). Both charges and vortices interact 
via short-range Yukawa potentials and behave thus freely when separated by distances 
larger then the scale (1/m). In presence of any dissipation mechanism (which would not 
alter the other two phases) this third phase corresponds thus to a metallic phase of the 
Josephson junction array 0. 

In conclusion we can represent the phase diagram of our model as follows: 

mlG{ml)7v ^ ^ [ I < 1 , superconductor {Za x Rb) , 



- > 1 , superinsulator (Ra x Zb] 
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( -g < mlGUl)n ' SUpcrCOuductor {Z A X i?^) , (4 23) 

'^^^^^>1^\^!^<I<'^^ metal (i?. X i?,) , 

l^ miGirniu ^ superiusulator (Ra x Zb) , 

where we have indicated in parenthesis the global symmetries of the various phases. In 
fig. ^ we plot the (numerically computed) function mlG{ml)n/ ^ for the value n = ln5. 
This gives an indication that a window for the metallic phase is open for ml just larger 
than 1, while in the regime ml < 0(1), relevant for Josephson junction arrays, a single 
phase transition from a superconductor to a superinsulator aX {e/g) = 1 is favoured. 

The experimental results for Josephson junction arrays are plotted in fig. ^ . These 
are essentially resistance measurements as a function of temperature in arrays with O(IO^) 
cells. The zero-temperature extrapolation of these results indicates a quantum phase 
transition between an insulator and a superconductor in the vicinity of the self-dual point 
Ej/Ec = 2/7^2 ~ 0.2. 



5. Quantum Hall phases and anyon superconductivity 



Recently it has been suggested that, in presence of Ug offset charge quanta per site 
and n^j, external magnetic fiux quanta per plaquette in specific ratios, Josephson junction 
arrays might have incompressible quantum fiuid fl^ phases corresponding to purely two- 
dimensional quantum Hall phases for either charges [1^] or vortices |T^ . 
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In analogy with the conventional quantum Hall setting one expects the charge 
and vortex transport properties to depend on the filling fractions {nq/nep) and (rifp/nq), 
respectively. Due to the periodicity of the charge-vortex coupling, however, (n^) is 
defined only modulo an integer as far as charge (vortex) transport properties are concerned. 
Using this freedom one can thus define effective filling fractions (we shall assume riq > 0, 
71^ > for simplicity): 

rtq 

= + , < i., < 1 , 

Us ^^--^^ 

= r 1- , r 1+ ' < < 1 , 

riq - [nq\ + [n^\ 

where [uq]^ indicate the smallest (greatest) integer greater (smaller) than Uq. These 
effective filling fractions are always smaller than 1. 



The masses of charges and vortices are given in (|3.19|) . For given values of Uq and 
which admit an incompressible quantum fiuid ground state one expects charges to bind 
vortices and form a charge quantum Hall phase in the regime where charges are heavier 
than vortices, i.e. e/g > 1. For given values of Uq and we thus expect a charge quantum 
Hall phase at filling Uq for e/g > 1 and a vortex quantum Hall phase at filling z/<^ for e/g < 1 . 
Correspondingly, these two regimes were analyzed in [l^] and [|ll|, respectively. In ||12|| , 
however, it was pointed out that e/g cannot be too small for the vortex quantum Hall 
phase, since for e/g '^1 the effective vortex band mass due to the periodic array becomes 
exponentially large and vortices loose their mobility. 

In the following we shall assume the existence of these quantum Hall phases and 
discuss them in the framework of the gauge theory representation of Josephson junction 
arrays in the self-dual approximation. The idea is as follows. For = = we have 
derived that the gauge theory describing Josephson junction arrays (in the self-dual ap- 
proximation) is given by ( |1 . 1| ) with periodic charge-vortex (mixed Chern-Simons) coupling 
and the identifications (|3.21|) . This gauge theory describes the dynamics of charge and 
vortex fiuctuations of the array in absence of external offset charges and fiuxes. Drawing 
on previous experience [|| with the quantum Hall effect we shall modify this gauge the- 
ory in order to describe charge and vortex fiuctuations about a homogeneous ground state 
with riq charges and vortices per plaquette, describing a quantum Hall fiuid for either 
charges or vortices. We shall then analyze how the periodicity of the new Chern-Simons 
charge-vortex couplings affects this picture. To this end we shall consider the Euclidean 
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partition function of the new gauge theories, enforcing the periodicity by appropriate topo- 
logical excitations and we shall study the resulting zero-temperature phase diagram. Given 
the expected jump in the relevant effective filling fraction at e/g = 1, we shall consider 
two different gauge theories in the regimes e/g > 1 and e/g < 1. 

5.1. Gauge theories for the quantum Hall phases 

Let us begin with the charge quantum Hall phase for e/ g > 1. To this end we consider 
the (Minkowski space-time) gauge theory with Lagrangian 

C, = -^F^F^ + -A^e^-^d^B, - ^f^F - ^F^^ + ^A^e^^^d^A, . (5.2) 

The main differences with respect to ( p..l| ) are the addition of a pure Chern-Simons term 
for the A^ gauge field and a new coupling term proportional to F^f^. The Gauss law 
constraint associated with the A^ gauge field now assigns a vorticity 

1 



2Z/gK2 



q (5.3) 



to a charge q = J dPx q^^'^^'^ (since all the gauge fields are massive there are no corrections 
to this equation from boundary terms). The F^f^ coupling then associates a corresponding 
magnetic moment fx oc fq/ ng"^ to these composites. We have also rescaled the coefficient of 
the mixed Chern-Simons coupling by a factor of 2 (compare with ( |1 . 1|) ) while maintaining 



the definitions (|3.21|) . This factor of 2 is a well-known aspect of Chern-Simons gauge 
theories [|2^ . Indeed, the vorticity ( ^.3|) has a back- reaction on the charges since it also 
couples to A^ via the pure Chern-Simons term. With our rescaling, the total current 
coupling to A^ is given by {2q'^^^'^^'^ -\- 2K^z^g^^°''*^^) and using ( |5T^ ) we see that the total 
"dressed" charge of the charge- vortex composite is indeed q. The rescaling ensures thus 
that dressed charges maintain their nominal value. 

The effective Lagrangian for the charge degrees of freedom, obtained by integrating 
out is given by 

>Cfff. = -T^ B^t^-'^d^B, + . . . , (5.4) 

where the ellipse stands for higher-derivative terms which are suppressed at long distances 
by inverse powers of a mass. Following Q] we introduce as external probes a conserved 
vortex current (j)^ and the electromagnetic gauge field A^™': 

3 

^fff. + ^r^^cW + = C^s. + ^ A'^^-e^-'^d^B, + K^B.r ■ (5.5) 

22 



Integrating also over the charge gauge field we find the effective Lagrangian 
C,s. [A-^-A^) = ^ A-^-e^^^d^A^-^- + K'l^q A-^-ct)^ + nt^i^^ ^p^e^^-^^p^ . (5.6) 

Prom this effective Lagrangian we learn two things. First of all, the electromagnetic re- 
sponse of the system is encoded in the induced current 

JL. = TT^^eff. {A--- A, = 0) = e^-d^Al--- , 

^' ' (5.7) 

where E is the applied electric field. This represents a Hall current with Hall conductivity 
given by 

^2 

an = —^q . (5.8) 

ZTT 

Secondly, the last two terms in ( ^.6|) tell us that represent charge and flixx carrying 
anyons with charge-fiux relation and fractional statistics given by 

q = UqK^4> , 

(5.9) 

9 = VqK^(j)'^ . 

An excitation carrying no effective vorticity can be obtained by combining a charge q with 
a vortex (f) = +q/2uqK'^ , so that the bare and induced vorticities cancel. This excitation has 
the standard electromagnetic coupling with q^ representing its conserved current. 

The Gauss law following from ( |5.6| ) then assigns to this excitation also a magnetic fiux 
qjvqK^ . All excitations in the model are therefore anyons satisfying (|5.9| ). Note that the 
magnetic moment can be written as /U oc 2S/Kg'^, where 5" = is the fractional spin 
associated with the fractional statistics (|5.y|). 

For Uq = p/n, with p and n coprime, n flux quanta 1/k have p units k of charge. Since 
in Josephson junction arrays the charge degrees of freedom are bosons (Cooper pairs), this 
excitation must also have bosonic statistics, i.e. 6 = even integer. This requires that pn 
must be an even integer. The allowed filling fractions are thus 

p 

q = — , pn = even integer , (5.10) 

n 



in accordance with 30 
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Note that there is no Chern-Simons term in the effective action for the vortices, 
obtained by integrating out S^. Indeed, the bare and induced Chern-Simons terms for 
Afj, cancel exactly. We thus conclude that (^3) is indeed the appropriate gauge theory to 
describe the charge quantum Hall phase of Josephson junction arrays. 

The gauge theory ( ^.2|) has a hidden duality, which can be made manifest by rewriting 
the Lagrangian as 



e' = 



(5.11) 



1 - 

g-^ q 



and defining a new gauge field B'^ = B^ + fqA^. Indeed, in terms of and S^, ( |5.11J ) 
coincides with ( |1 . 1| ) upon substituting e e' (and k 2k). In the sector in which 
f^J, + '^qF^ = the only kinetic term of ( [5.11| ) is contained in — (l/2e'^)-F^F^. Therefore 
niqij) = l//^e'^ is the mass of the anyonic charge-fiux composites. The gap for collective 
oscillations is given by the modified topological Chern-Simons mass 

M, =m(e',9,2K =-2- = ^ . 5.12) 

In the representation ( p.llj ) it is also manifest that our modified gauge theory can be 
defined only in the range 

!<-< — . (5.13) 

9 ^q 

For e/g l/z^^ the anyon mass rrig^ vanishes, while the topological mass Mg diverges. 
The gauge theory describing the vortex quantum Hall phase is the dual of (|5.2|), 



= -T^F^F^ + -A^e^^'^d^B, - ^f^F - ^F^r + ^B^e^^'^d^B, , (5.14) 
2e^ n Ig^ n 

and contains a pure Chern-Simons term for the gauge field B^ and the corresponding 
magnetic moment interaction. Again, the rescaling of the mixed Chern-Simons coupling 
by a factor of 2 ensures that dressed vorticity maintains its nominal value. 

In this case there is no pure Chern-Simons term in the effective action for the charges, 
while the vortex effective Lagrangian is given by 

^ts = -j^A^e^^^d^A, + . . . , (5.15) 
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at long distances. We probe the vortex response by coupling (|5.15|) to a gauge field G^^^- 
such that e^^'^QaCf 

= Qext. describes an external current distribution. We also introduce 
an additional conserved charge current coupling to in order to probe the quantum 
numbers of excitations: 

^ts + Xcf-e'^'^-d^A^ + -^A^q'^ . (5.16) 
Integrating over the vortex gauge field we find the effective Lagrangian 

AfT. {Gf-,q,) = ^ Gf-e^-'^d^G^- + ^ G^^'-q^ + ^ q,e^'"'^q. • (5.17) 

This leads to the following induced vortex current: 

- ^ ^^^eff. (G--, = 0) = ^ Ct. ■ (5.18) 

This equation embodies the quantum Hall effect for vortices. While charges react to 
external electromagnetic fields, vortices react to external electric currents: in the normal 
case the induced vortex current is perpendicular to the applied electric current, in the 
quantum Hall phase it is parallel, with coefficient proportional to z^^. From the last two 
terms in (|5.17| ) we read off the fiux-charge relation and the fractional statistics of the anyon 
excitations represented by q'^: 

,2 ' 



(5.19) 



^2 ■ 

For Vff) = p/n, n charge quanta k carry p quanta 1/k of fiux. Since vortices are also bosons, 
we find the same allowed values of z^<^ as for the charge quantum Hall phase: 

p 

= — , pn — even integer . (5.20) 

n 

The self-dual representation analogous to ( p.ll|) is given by 

% ^ (5.21) 

9 ~ 

and coincides with ( |1.1| ) upon introducing a new gauge field A'^ = A^ + z/^-B^ and sub- 
stituting g ^ g' (and k 2k). In this case m^g = 1/l^g'^ is the mass of anyonic 
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flux-charge composites, while the gap for collective oscillations is given by the topological 
Cher n- Simons mass 

^ m(e, g\ 2k) = = . (5.22) 



Clearly, ( ^.21|) is defined only in the range 



e 



z^w, < - < 1 . (5.23) 
9 

Again, for e/(7 ^ z^^ the anyon mass m^g vanishes, while the topological mass diverges. 
Combining (|5.13[ ) and ( p.23| ) we find the overall condition 



u^, < e/g < , (5.24) 

which we interpret as the regime in which a homogeneous ground state with Uq charges and 
Ufj) vortices per plaquette can exist. Presumably, for e/^r < z/<^ and e/g > l/i^q the ground 
state consists of an Abrikosov-type cristal for charge-flux composites. In particular, ( |5.24| ) 
tells us that in Josephson junction arrays Ec / Ej cannot be either too large or too small 
for the existence of quantum Hall phases. Although its origin is different, this condition 
agrees with the result of []12| (at least for the vortex quantum Hall phase). 



5.2. Periodic Chern-Simons terms and phase structure analysis 

In the following we shall analyze how the above picture is modified when we impose the 
distinctive feature of Josephson junction arrays, namely the periodicity of charge-vortex 
couplings, encoded in our formalism in the Chern-Simons terms. This is achieved by 
introducing appropriate topological excitations in the Euclidean lattice partition functions 
of (ISHD and iWM )- 



Let us begin with the gauge theory (|5.11| ) for the charge quantum Hall phase. Its 
Euclidean lattice partition function coincides with (|3.20| ) upon substituting e — > e', — > 



S-^ = -f VqA^ and rescaling the mixed Chern-Simons term by a factor of 2. Therefore 
we present here only the coupling of the topological excitations enforcing the periodicity 
of the mixed Chern-Simons term A^k^j^B^ and the Wilson and 't Hooft loops (|4.6| ): 



Sq = S^... + ilp^A^ {Q^ + M^) + iln^B^M^ -f il\A^q^ -f il(j)K^ B^(l)^ , (5.25) 

K2 
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where we have used the representation fq = p/n. Due to the change — > the 
periodicities of the two original gauge fields are changed from ( p.25|) to 



B^^ B^ + ^b^ , b^ez 



(5.26) 



The displayed terms in ( p.25|) can be rearranged as follows: 

Sq = J2--- + ilpV^A^ (q^ + -^q^ - '^4^^ + iln^Bl (m^ + ^</)^') , (5.27) 

so that the whole model is reformulated in terms of the gauge fields and and we 
can use thus the results of the previous section. In particular we obtain 



with 



Top 



X 1 1 

+ i— — M 

^ I ^'^V2(M|-V2)^' 

and Mq defined in (|]T|). 

At this point we can repeat verbatim the analysis of section 4. The phase structure of 
( ^.11|) with periodic Chern-Simons term is governed by the topological excitations and 
M^. The phase in which both these topological excitations are dilute corresponds to the 
charge quantum Hall phase discussed above, which constitutes an incompressible fluid of 
charge-fiux composites with short-range interactions. The stability of this phase depends 
entirely on the condensation conditions for the two types of topological excitations. If 
condenses we obtain a phase in which charges q = up (and multiples thereof) are completely 
screened, while fiuxes (p are logarithmically confined: this is a conventional superconducting 
phase with a charge np condensate. Using the trick (|4.1UD we can identify the Coulomb 



law for fiuxes with the London form of the electromagnetic response. If condenses 
we obtain, instead, a phase in which excitations with quantum numbers {q/{pK) — (pn/n) 
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interact logarithmically. This means that the only non-confined excitations in the model 
must carry both charge and fiux in the combination 

= =^ — = K fq ■ (5.30) 

np n (p 

These excitations are completely screened, while all other combinations of quantum num- 



bers are logarithmically confined. This (logarithmic) oblique confinement [|T6[ phase 



describes thus a charge-flux superfluid phase. Since the condensed composites carry charge, 



this is actually an any on superconductivity phase ||17|| with a charge pn and vorticity u/k 
condensate. Indeed, the electromagnetic response, obtained again by the trick (|4.10|) has 
still the London form. 

We now have to study the range of parameters in which these three phases are realized. 
The analysis goes exactly as in section 4, giving the result 



< 1 - ^ e' 1 



^ < — , conventional superconductor 
> — , anyon superconductor , 



J < ^ ^ , conventional superconductor , 

X„ > 1 — > < < — < ^ , charge quantum Hall phase , (5.31) 



'■9 ) l^qXq g 

^ > , anyon superconductor , 

^ _ MqlG (Mql) TT pH 

p 2 

It is now harder to disentangle the phase diagram in terms of the original (Josephson 
junction) parameters e and g since e', and consequently Mg depend themselves on the 
ratio e/g. 



As was pointed out in [|T2|, the periodicity of charge- vortex couplings is the distinctive 
feature of Josephson junction arrays which allows effective filling fractions of order 0(1) 
and which is therefore expected to favour the formation of charge and vortex quantum 
Hall phases. From the above result, however, it is clear that the same mechanism can also 
destabilize these quantum Hall phases as follows. The condition for a charge quantum Hall 
phase is given by Xq > 1. The filling fraction parameters p and n enter this condition in 
two ways. First there is the explicit dependence of Xq on the product pn; secondly there is 
the dependence of the gap Mq on the ratio Uq = p/n. The former favours filling fractions 
with a large product pn; too large numerators and denominators are however presumably 
suppressed by the same mechanism as in the conventional quantum Hall effect. The latter 
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has the foUowing effect. Decreasing both Uq and e/g makes both the gap Mg and the 
ratio e' / g smaUer; the parameter Xq approaches a fixed value, while e'vqj g can decrease 
indefinitely till it becomes favourable for the system to expel the magnetic fiux and form 
a conventional superconductor. Increasing both Vq and ej g, instead, makes both the gap 
and the ratio e' j g larger; for large values of M^/ the parameter Xq tends to zero (see fig. |l]), 
while e'vql g can grow indefinitely and it quickly becomes favourable for the charge-fiux 
fiuid to condense into a superfiuid, so that the system becomes an anyon superconductor. 
For values evqj g > 1 we would expect a charge-fiux cristal. 

The analysis of the gauge theory (|5.21| ) for the vortex quantum Hall phase follows ex- 
actly the same steps with all "electric quantities" and "magnetic quantities" interchanged. 
Therefore we present here only the final result: 



< 1 



< Uff) , anyon superconductor 
> U(f) , superinsulator , 



fr < 1 anyon superconductor , 

Xfj, > 1 ^ ^ ^ P ^ , vortex quantum Hall phase , (5.32) 
> UfjjXfj) , superinsulator , 

M^IG{M^I)tx pn 



Xrh = 



where is defined in ( |5.22| ) and we have used the representation z^^ = p/n. Starting from 
the vortex quantum Hall phase and decreasing both g/e and makes both the gap 
and the ratio g' /e smaller: it becomes eventually favourable for the system to expel the 
offset charges and become a superinsulator. In real Josephson junction arrays this phase 
would presumably be an insulating Abrikosov-type cristal of charges due to the small but 
finite ground capacitance Cq. Increasing both g/e and z/^ makes both the gap and the 
ratio g'/e larger; the quantity l/X^f, tends to infinity (see fig. |l|) while e/g' decreases and 
the fiux-charge fiuid of the vortex quantum Hall phase condenses again into a fiux-charge 
superfiuid, becoming thus an anyon superconductor. For even smaller values of e/g < 
we would expect again a fiux-charge cristal. 



6. Three dimensions 

In this section we generalize our results (for zero offset charges and external mag- 
netic fiuxes) to three-dimensional Josephson junction arrays. While these are not (yet) 
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experimentally accessible, we find it nonetheless interesting to construct their gauge the- 
ory representation and to study the differences with the two-dimensional case. Clearly, 
self-duality is lost in three dimensions since the fluctuating degrees of freedom are charges 
and closed vortex loops. However, it is still possible to introduce the approximation which 
allows a coupled gauge theory representation, as we now show. 

Up to eq. (|3.10D the analysis parallels exactly the two-dimensional case. The decom- 



position analogous to (|3.10| ), however, requires the introduction of the three-dimensional 
generalizations of the lattice operators K^^, and K^^. These are given by the three-index 
lattice operators 



(6.1) 



where Sa and 5"^ are the shift operators (2.9) and and An are the finite difference 



operators ( |2.15| ) . As in the two-dimensional case, these two operators are interchanged (no 
minus sign) upon summation by parts on the lattice. Moreover they are gauge invariant, 
in the sense that they obey the following equations: 



pAjy K pA p ApK^jyp , 

flU pA p ApKpjjp AjyK pi, p . 



(6.2) 



Finally they satisfy also the equations 

K^^pKpXu = - {Spx5,y^ - 5pujSv\) A + (^5p\A„Aaj - 5„xA^A^ 

+ (^5„^ApAx - 5p^A„Ax) , (6.3) 
piypK pi/ij K pi/ pK pi/i_j 2 ^6pi_jA ApA^^ 

Using these operators we can decompose Vp as 

Vp = ApUi + Apa + Kpo,pipai3 , (6.4) 

with m E Z and \a\ < 1. The V'a/j's are restricted by the fact that the antisymmetric 
combinations qp^ = Kp^aVa = Kpi^aKaXpi/^xp must be integers. We can thus trade the 
original sum over the four independent integers {vp} for a sum over the seven integers 
{m,qp,y} subject to the constraint Apqp^ = A^qpi, = 0. This constraint eliminates the 
three longitudinal degrees of freedom of qp^, so that {m,qp,^} with the above constraint 
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describes only four independent integers. After shifting the $ integration domain using 
the sum over {m} and performing the resulting trivial $ integration we are left with 

(6.5) 



2 
PI 



21qEj 



The constraints are solved by introducing a real antisymmetric gauge field b^^, and 
an integer gauge field a^: 

Pfj. = K^apbaP , baP G R , 

(6.6) 

Repeating the same steps as in the two-dimensional case we find 

^^^^ ^ (6.7) 

^ 1 J)- 

^ -A IIqUij 

which is the three-dimensional analogue of ( |3.13| ). Here, K^^pbctp maintains its interpre- 
tation of the conserved four-current of charge fluctuations, while K^^aO-a. represents the 
fluctuations of closed vortex loops. 

Since the magnetic fluctuations are represented by closed vortex loops we cannot 
render the partition function self-dual, as in the two-dimensional case. However, it is still 
possible to introduce a bare kinetic term for the vortex loops with a coefficient tuned so 

2 

that it can be combined with the Coulomb term for the charges into 2NH e 'iij'iij- 

As 

in two dimensions we have to introduce new integers via the Poisson summation formula to 
guarantee that the charge current K^ctpbap remains an integer. In three dimensions these 
are two-index antisymmetric integers M^jj satisfying the constraint A^M^j, = K^M^j^ = 0: 



Z= Yl / ^^^M exp(-^) 



{Q^} 
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After a rescaling 





27r 






Boi — ■ 




J — Oni 


VnIoI 



A - 2^ 

VNl 



(6.9) 



we obtain finally 



{Qm} 

Q Y-W^ P ^ .W^^, , R Jol' (6-10) 

X ^ 

+ {loQoAo + lQ^Ai) + {IqIMq.Bq^ + IoIMmBm + l^M.jBij) , 

where k^^p and k^^p are the analogues of defined in terms of derivative operators 
rather than finite difference operators (and satisfying the correspondingly modified eq. 
(|673|)) , Fij and fi are the spatial components of 

F — h A 

_ 1 (6.11) 

and the coupling constants (dimensionless) and g'^ (with dimensions mass^) are given 
by 



e^ = 2NlEc, n = N, g^ = l-Ej. (6.12) 
The plasma frequency is given again by a product of these coupling constants, 

m = ^ = ^/2N^EcEj , (6.13) 

TT 

and for ml < 0(1) we can choose Iq = I. Also in three dimensions we have thus obtained 
a coupled gauge theory in the limit of infinite magnetic permeabilities, encoded in the 
absence of the time components of both Ff^^, and in the kinetic terms. As in two 
dimensions we shall henceforth consider the relativistic limit of this gauge theory: 



2 



Z= J2 j^^^ j exp(-5) 



{Qm} 



(6.14) 
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This is a pure gauge theory representation of a model of interacting charges and closed 
vortex loops. The identification of the physical degrees of freedom is analogous to the 
two-dimensional case: 



3 

K2 



charge — _ U „R „ 

J vortex _ ^ t A 
Vfiv — - 1 "^MfQ^c 



(6.15) 



27rK2 

The model (|6.14| ) is a (Euclidean) lattice version of the so called BF gauge theory [|T^ 



whose Lagrangian is given by 

Ik 1 

/^.p = d^B^p + d^Bp^ + dpB^^ , (6.16) 

Here describes an ordinary (3+l)-dimensional photon with field strength given by 
F^^ = d^A^ - d^Af, and dual field strength F^"' = ^e^^^^PFxp. This photon has a 
topological BF coupling to an antisymmetric Kalb-Ramond gauge field -B^j^, whose 



field strength is given by the three- form -F^jyp. The first term in ( 6.16 ) represents the 
kinetic term for the Kalb-Ramond gauge field. In addition to the usual invariance un- 
der gauge transformations of (|6.16| ) is also invariant under gauge transformations 
B^i, — » B^i, + dp,Kj, — d^A^. The Kalb-Ramond gauge theory describes a single mass- 
less scalar degree of freedom since all the components Bio are Lagrange multipliers and 
Bij has only one transverse component. When coupled to the usual Maxwell theory via a 
BF term, the Kalb-Ramond sector induces a topological mass 

m = (o.lTj 

for the photon. As already pointed out in ( |6.13|) this mass represents the plasma frequency 
of Josephson junction arrays. The BF system is thus the natural three-dimensional gen- 
eralization of (|LT|), as expected. 

As in two dimensions, the integer-valued variables Qp and Mp^, appearing in ( |6.14| ) 
represent (Euclidean) topological excitations whose role is to make the charge-vortex BF 
coupling periodic. They satisfy the constraints 

dpQp = , 

. . (6.18) 

dpMp, = d^Mp, = . 
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The "electric" topological excitations are exactly as in two dimensions; the "magnetic" 
topological excitations, instead, describe compact surfaces on the lattice (in a dilute phase) 
or infinite surfaces (in a dense phase). 

The phase structure of three-dimensional Josephson junction arrays (in our approx- 
imation) is thus determined by the statistical mechanics of a coupled gas of lattice loops 
and surfaces. Its partition function can be easily obtained by a Gaussian integration over 
the gauge fields and B^j, in ( |6.14| ): 

^Top = ^ exp (-5'top) , 

St,, = E W Q.'J^.Q- + 4^ Af„„^i4|^i|f^M„, (6.19) 



X 



2/2 ^'^^2_v2^'^ '2 m2-V2 



with m defined in (|6.17| ) . This partition function can be interpreted as the Euclidean parti- 
tion function for a lattice model of interacting particles (whose world-lines are parametrized 
by the closed loops Q^) and closed Nielsen- Olesen type strings (whose world-sheets are 
parametrized by the compact surfaces M^^). In a derivative expansion the string action 
takes the form 

^strings = Yl —Mf^uM^. + ... . (6.20) 

f\jC 

X 

In the dilute gas approximation, where M^i, can take only the values 0, ±1, this term 
measures the area of the world-sheet and is thus the standard Nambu-Goto term ||^, with 
string tension tt^/kc^/^ (remember that was the vortex mass in two dimensions). 

The parameter l/l plays thus the role of the Higgs mass in our lattice model; higher order 
terms in ( |6.20| ) involve both the curvature and internal excitations of the string. Particle- 
string interactions are encoded in the topological Aharonov-Bohm term (third term in 
( |6.19|) ) measuring the linking of the closed world-lines of particles and compact world- 
sheets of strings in four Euclidean dimensions. As in the two-dimensional case this term 
vanishes for loops and surfaces separated on distances much larger than 1/m. In this case 
the denominator reduces to m^V^ and, by using either one of the representations 

_ . (6.21) 
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and the equations ( |6.3|) one recognizes that the whole term reduces to (i27rinteger), which 
is equivalent to 0. 



Unfortunately, the statistical mechanics of random surfaces [^T| is much less under 



stood than its random loop counterpart and we cannot use self-duality arguments 



anymore. In order to proceed further we shall assume that the same three types of phases 
as in the two-dimensional case can exists and we will point out the differences that can 
nonetheless arise. First of all we can repeat the same argument as in the two-dimensional 
case to find that in the electric condensation phase the global gauge symmetry associated 
with Afj^ is spontaneously broken Ra Za while in the magnetic condensation phase it 
is the Kalb-Ramond global gauge symmetry which is spontaneously broken Rb — > Zb- 
Secondly we can consider the expectation values of the Wilson loop ( [4.6|) and the 't Hooft 
surface 

Sh = exp , (6.22) 

where (f)^,^ vanishes everywhere but on the plaquettes of a compact surface, where it takes 
the value 1. These expectation values are given by 

^^^^ ^ ^Top (Q^ + ^g^,M^.) 

ZTop{Qfi,M^„) ' (6 23) 

^g^s^ ^ Ztop (Qm^ ^t^'^ + ^"^(l^fiu) 

Ztop {QfjL, M^j,) 

With exactly the same computation as in the two-dimensional case we find the dominant 
contributions 

2kP ^^-V 



{^w)nia,g. cond. = exp [ ^ qfi,—=^qu 



o . 9 X (6-24) 



{SH)e\. cond. = exp I (pfi^— 

Small loop (surface) corrections can alter only the coefficients of the Coulomb potentials 
in these results. The long-range nature of the interaction kernels is associated with the 
Goldstone bosons due to the spontaneous symmetry breaking present in both the electric 
and the magnetic condensation phases. In four Euclidean dimensions the results ( |6.24| ) 
for the 't Hooft surface in the electric condensation phase implies that the self-energy 
of a circular vortex loop of radius R is proportional to RlnR. As in two dimensions 
this is tantamount to logarithmic confinement of magnetic fluxes and we conclude thus 
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that the electric condensation phase is stiU a superconducting phase. The result ( |6.24| ) 
for the Wilson loop in the magnetic condensation phase, instead, represents a perimeter 
law, implying a 1/r Coulomb potential between charges. The amount of energy required 
to separate a charge-anticharge pair is finite, although the interaction potential is long- 
range. We identify this as an insulating phase (as opposed to the superinsulator in two 
dimensions). Clearly, the dilute phase for both topological excitations corresponds again 
to a metallic phase. As already mentioned, it is harder to estimate the position of the 
phase transitions in three dimensions due to the lack of self-duality. 
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Figure Captions 



Fig. 1. The function "^''^(^"^^)^ for ^ = ln5. 

Fig. 2. Pliase diagram of fabricated Josephson junction arrays (adapted from the last 
paper in 0). Solid squares denote a transition from metallic behaviour to su- 
perconducting behaviour when the temperature is lowered; open squares denote 
a corresponding transition from metallic to insulating behaviour. 
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